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Based on the Hugenholtz-Van Hove theorem, it is shown that both the symmetry energy Esj,m(p) 
and its density slope L{p) at normal density po are completely determined by the global nucleon 
optical potentials that can be extracted directly from nucleon-nucleus scatterings, (p,n) charge ex- 
change reactions and single-particle energy levels of bound states. Adopting a value of m* /m = 0.7 
for the nucleon effective k-mass in symmetric nuclear matter at po and averaging all phenomeno- 
logical isovector nucleon potentials constrained by world data available in the literature since 1969, 
the best estimates of Esym{po) ~ 31.3 MeV and L{po) — 52.7 MeV are simultaneously obtained. 
Uncertainties involved in the estimates are discussed. 
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Nuclear symmetry energy Esym (p) , which encodes the 
energy related to neutron-proton asymmetry in the equa- 
tion of state of nuclear matter, is a fundamental quan- 
tity currently under intense investigation in both nuclear 
physics and astrophysics [l-16|. Both the magnitude 
and density dependence of Egymip) are critical for un- 
derstanding not only the structure of rare isotopes and 
the reaction mechanism of heavy-ion collisions, but also 
many interesting issues in astrophysics, such as the struc- 
ture and composition of neutron stars. Despite much 
effort made both experimentally and theoretically, our 
current knowledge about the Egymip) is still rather poor 
even around the saturation density of nuclear matter pg. 
More specifically, near po the symmetry energy can be 
characterized by using the value of Esymipo) and the 
slope parameter i(po) — Spo- ^'^'""'*''^'' 



of neutron stars, much more effort has been devoted re- 
cently to extracting the value of L from studying vari- 
ous phenomena/observables in terrestrial nuclear labora- 
tory experiments. These include the isospin diffusion, 
neutron/proton ratio of pre-equilibrium nucleon emis- 
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The Egymipo) is known to be around 28 — 34 MeV mainly 
from analyzing nuclear masses within liquid-drop models 



sions and isoscaling in heavy- ion reactions [111 . |33| - |37| . 
energy shift of isobaric analogue states [38| , neutron-skin 
of heavy nuclei 24 , 39| . and giant dipole as well as pygmy 
dipole resonances [40|-|42| . Unfortunately, the values of L 
extracted so far from these studies scatter between about 
20 to 115 MeV. Since all these phenomena/observables 
are in someway at least indirectly and qualitatively re- 
lated to the L parameter, it is very useful to know if one 
can directly express the L in terms of some relevant parts 
of the commonly used underlying nuclear effective inter- 
action. In this Letter, based on the Hngenholtz-Van Hove 
(HVH) theorem [i^l, it is firstly shown analytically that 
the i(p) is determined completely by the momentum- 
dependent global nucleon optical potential that can be 
directly extracted from nucleon-nucleus and (p, n) charge 
exchange reactions. By averaging all phenomenologi- 
cal global isovector optical potentials constrained by the 
world data available in the literature since 1969, the best 



[ITHia. The exact value of Esym{po) extracted in this estimates of E,y„(po) = 31.3 MeV and i(po) = 52.7 



way depends largely on the size and accuracy of nuclear 
mass data base as well as whether/what surface symme- 
try energy is used in the analysis. The empirical value of 
Esymipo) has been used in calibrating many-body the- 
ories with various interactions. The density slope i(po) 
has also been extensively studied but remains much more 
uncertain. Since its exact value is particularly impor- 
tant for determining several critical quantities, such as, 
the size of neutron-skin in heavy nuclei |20l425|. loca- 
tion of neutron dripline [i^, core-crust transition den- 
sity 0,13, H 27 -31 1 and gravitational binding energy 32 1 



MeV are obtained simultaneo usly . 
The general HVH theorem 
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is a fundamental relationship among the Fermi energy 
Ep, energy density ^ = pEav where Eav is the average 
energy per nucleon and the pressure of the system P at 
zero temperature. The HVH theorem has been strictly 
proven to be valid for any interacting self-bound infinite 
Fermi system and is independent of the precise nature 
of the interaction used. At normal density po where the 
pressure vanishes {P = 0), the general HVH theorem 
reduces to the special relation Ep ~ Ea^. According to 
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the general HVH theorem, the Fermi energies of neutrons 
and protons in nuclear matter of isospin asymmetry S = 
(Pn ^ Pp)/p an arbitrary density p are, respectively. 
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where t{k^^^) and C/„/p are the neutron/proton kinetic 
energy and single-particle potential, respectively. The 
Fermi momenta of neutrons and protons arc kp ~ 
kpil + Sy/^ and k'^p = kpil - 5)^/^, respectively, with 
kp = (3t:^ p/2)y^ . According to the well-known Lane 
relationship [4^, the Un/p can be well approximated by 
U„/p{p,5,k) ~ Uo{p,k)±Usym{p,k)S, where Uo{p,k) and 
Usym{p,k) are, respectively, the isoscalar and isovector 
(symmetry) nucleon potentials. While the Usym{p,k) 
still has significant uncertainties the Uo{p,k) especially 
at po has been relatively well determined P, 0, l^- 
By subtracting the Eq.Q from Eq.® and then mak- 
ing a Taylor's expansion in 5 of both the t{k^^^) and 
i/p 



Un/p{p,5, kp ), one can show that [45l. l46l| 

Esymip) = \ ^^^QI^°^ \kFkF + ^UsyjniP^kp). (5) 

Moreover, by firstly adding the Eq.© with Eq.(|3]) 
and then expanding in S again both the t{kp^^) and 



left hand side relates to the momentum slope of single 
nucleon potentials when both the t{kp^^) and Un/p are 
expanded as power series of 6 through the neuron/proton 
Fermi momentum. 

To evaluate the symmetry energy and its density slope 
at pq using the Eqs.((5|) and we need to know 

the momentum dependence of both the J7o(poj^) ^^nd 
UsymipoT k). This information can be obtained from the 
energy dependence of the depth of the real part of the 
Global Optical Potential (GOP). The latter is normally 
parameterized in the Woods-Saxon form. Over the last 
few decades, great progress has been made in developing 
the unified GOP for both nuclear structure and reac- 
tion studies. At negative energies the GOP can be con- 
strained by single-particle energies of bound states while 
at positive energies it is constrained by nuclear reaction 
data [irl - figt . In the present study, we shall take the full 
advantage of the systematics accumulated in this field. 
The momentum dependence of Uo{p, k) is conventionally 
taken into account by using the nucleon effective mass 
approximation. In terms of the effective k-mass m* de- 



fined by m*/m = [1 + 
L{p) can be rewritten as 
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Un/p{p, S, k^p) while noticing the Eq.®, we obtain for For the present study, we only need to know the m* and 

information about the Usym{p,k) at po and the corre- 
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It is worth mentioning that the analytical expressions 
in Eq.(IS]) and Eq.® are valid at any density and the 
values of both Esy,n{p) and L[p) can be easily calcu- 
lated simultaneously once the single-particle potential is 
known. The most critical advantage of the expression 
in Eq.® is that it allows us to determine the L{p) di- 
rectly from the value and momentum dependence of the 
single nucleon potential at p. Essentially, this enables 
us to translate the task of determining the density slope 
of the symmetry energy into a problem of finding the 
momentum dependence of the Uo{p,k) and Usym{p,k). 
Very fortunately, the latter is accessible at least at nor- 
mal density from phenomenological optical model analy- 
ses of nucleon- nucleus scatterings, (p, n) charge exchange 
reactions and the single-particle energy levels probed by 
pickup and stripping reactions [47H49[. The fundamen- 
tal reason for this interesting feature can be understood 
mathematically from inspecting the Eqs.([3]) and (j4]). It 
is seen that the right hand side of these equations relates 
to the density slope of the symmetry energy, while the 
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spending Fermi momentum {kp = 1.36 fm~^). For the 
nucleon effective mass, we adopt the value of m* /m = 
0.7 ± 0.05 [H^ widely used in the literature, see, e.g. 
[sH . This value is consistent with the Uoipo,E) = 
— (50.0 — 0.30i?) from the global optical model analysis of 
nucleon-nucleus scattering experiments 14 7|. 1521 and the 
microscopic many-body calculations |53l. |54|| . ^Since the 
Uq from data analyses is usually expressed as a function 
of nucleon energy E, a dispersion relation determined by 
Uo{po,k) = -{50.0 - 0.30[t(fc) + f/o(po,fc)]} has to be 
used to obtain the m* /m ~ 0.7. 

The symmetry potential UsymiPoi E) can be deduced 
from studying the i sosp in dependence of the optical po- 
tential 43, |48|, |55|, |56| using 1) elastic scattering of a 
proton and a neutron from the same target at about 
the same beam energy after correcting the Coulomb en- 
ergy for protons; 2) proton scattering at the same en- 
ergy on several isotopes of the same element; 3) (p, n) 
charge exchange reaction between isobaric analog states. 
Since the 1960s, there are many sets of GOPs deduced 
from phenomenological model analyses of the world data 
existing at the time [H^, H^l- While some of the anal- 
yses assumed an energy independent symmetry poten- 



tial, see, e.g., j47l |57| . a significant number of stud- 
ies considered the energy dependence [HI]. We also 
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well with that deduced from the experimental binding en- 
ergies of nuclei. Similarly, a value of L{pq) = 52.7 MeV 
is obtained from Eq.®. 
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FIG. 1: (Color online) Energy dependence of the nuclear sym- 
metry potential Usym{po, E) at normal density from different 
global optical model fits. The solid lines are in the energy 
ranges where the original analyses were made while the dashed 
parts are extrapolations. 
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notice here that there are also some local/global phe- 
nomenological optical potentials with mass-dependence 
and semi-microscopic optical potentials involving com- 
plex nuclear structure calculations (see Ref.jH^ and ref- 
erences therein). Here we only use results from energy 
dependent studies. In these analyses the symmetry po- 
tentials are normally described by using a linear function 



of the form Usymipo,E) = Gsym 

1 are all the energy dependent symmetry potentials that 
we can find in the literature |58| . While some of them 
were extracted from data in the same energy range, some- 
what different kinds of data and models were used. All 
together, they represent parameterizations constrained 
by the world data on nuclcon- nucleus scatterings, (p, n) 
charge exchange reactions and single-particle energies of 
bound states. Assuming that these various global energy 
dependent symmetry potentials are equally accurate and 
all have the same predicting power beyond the original 
energy ranges in which they were analyzed, an averaged 
symmetry potential of 



FIG. 2: (Color online) The magnitude of each term in the 
nuclear symmetry energy Esym{po) and its density slope L{po) 
at normal nuclear density. 

At the mean-field level, isospin effects are determined 
by the quantity Usym ' ^- Since the latter is generally 
much smaller than Uq [ssI . [60j . it is no surprise that 
there are large discrepancies among results from differ- 
ent analyses. There are thus also probably large error 
bsyrriE. Shown in Fig. bars associated with our best estimates for the Egymipo) 

and L{po). However, due to the wide range of experi- 
ments involved and the various methods used in analyz- 
ing the data, it is impossible for us to calculate accu- 
rately these error bars. Nevertheless, it is instructive to 
estimate the uncertain ranges of each term in calculating 



the Esymipo) Si,nd L(pq). By using the ranges of the pa- 
rameters asym and bsym shown in Fig. 1 , the uncertainties 



of all terms contributing to the Esym{po) 
obtained and shown in Fig. 2. The Es. 



and L{po) are 
.(1) - 1"''^ 



Usym{po,E) = 22.75 -0.21E 



(9) 



is obtained (thick solid line in Fig. 1). Besides deter- 
mining the values of Esym{po) and L{po), we expect this 
expression to be useful for calibrating Usym (p, E) at ab- 
normal densities predicted by various many-body theo- 
ries. Under the assumptions mentioned above, it repre- 
sents the best fit to the global symmetry potentials con- 
strained by the world data up to date. With this best 
estimate for the energy dependent symmetry potential at 
po, it is then straightforward to get the best estimates 
for the Esym{po) and L{po) using Eq.([7]) and Eq.®, 
respectively. Because E = ~\Q MeV at po, the cor- 
responding symmetry potential is then Usymipo^kp) = 
22.75 - 0.21 X (-16) = 26.11 MeV. With this we obtain 
the Esymipo) ~ 31.3 MeV. This value agrees remarkably 



^symy.-^/ ~ 3 2„i' 

denotes the kinetic energy term with the effective mass 
TO* and the Esym{'^) = \Usym{po,kF) is the symmetry 
potential contribution. It is seen from Fig. 2 that the 
two terms are comparable. The uncertainty for calcu- 
lating Egyra{^) is due to the effective mass to*. With 
the m* /m = 0.70 ± 0.05 we adopted [H^, an error bar 
of ±1.0 MeV is obtained for the i?si/in(l), as marked 
by the red box. From the results shown in Fig.l, the 
symmetry potential Usymipo^kp) at the Fermi momen- 
tum kp is found to lie in the range of 26.11 ± 7.0 MeV. 
So, put all together the error bar of Egymipo) is ap- 
proximately ±(1.0 + 3.5) = ±4.5 MeV. The L(po) has 

three terms L{1) = l^:^, L{2) = '^UsyraiPQ-.kp) and 
i(3) = Q^'«^^(P'fe) |fe^fcj,. While both the L{1) and L{2) 
are positive, the L(3) is negative because of the decreas- 
ing symmetry potential with increasing energy. Similar 
to estimating the error bars for the Esym{l) and Esym{2), 
an error bar of ±2.0 MeV is estimated for L{1) and ±10.5 
MeV for L{2). To evaluate the error bar for L(3), we es- 



4 



timated an uncertainty of 0.21 ± 0.1 for &sym from Fig.l. 
Then, an error bar of ±10 MeV is obtained for L(3). 
Thus, the total error bar of L{pq) is approximately ± 
22.5 MeV. 

In summary, using the Hugenholtz-Van Hove theorem, 
it is shown that both the symmetry energy and its den- 
sity slope at normal density are completely determined by 
the global nucleon optical potentials. Using the nucleon 
effective mass m*/m = 0.7 and the nuclear symmetry 
potential Usym{po, E) = 22.75 — 0.21^^ obtained from av- 
eraging all phenomcnological isovector nucleon potentials 
constrained by world data up to date, the best estimates 
of the symmetry energy and its density slope at po are 
found to be Esym{po) = 31.3 MeV and L{po) = 52.7 
MeV, respectively. Uncertainties involved in the esti- 
mates are discussed. 
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